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Abstract 

We describe a general search for resonances decaying to a neutral e/i final state in pp collisions 
at a center-of-mass energy of 1.96 TeV. Using a data sample representing 344 pb"-^ of integrated 
luminosity recorded by the CDF II experiment, we compare Standard Model predictions with the 
number of observed events for invariant masses between 50 and 800 GeV/c^. Finding no significant 
excess ( 5 events observed vs. 7.7ib0.8 expected for Mg^ > 100 GeV/c^ ), we set limits on sneutrino 
and Z' masses as functions of lepton family number violating couplings. 

PACS numbers: 11.30Fs,12.60Cn,12.60Jv,13.85Rni 
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An observed excess of high- mass opposite sign electron- muon pairs at the Tevatron would 
provide evidence of physics beyond the Standard Model (SM). Nearly every extension of the 
SM predicts flavor changing neutral current effects. General supersymmetric models, for 
instance, permit the violation of R-parity symmetry (RPV) and describe the lepton family 
number violating (LFV) production and decay of sneutrinos (u), the scalar superpartners 
of SM neutrinos Models with additional gauge symmetry can also accommodate an 
e/i signature through LFV decays of a new heavy neutral gauge boson, the Z' ^. Signals 
from such processes may be easily detected at the Tevatron. SM backgrounds are small 
and characterized by invariant mass spectra that lie well below the range presumed for new 
physics. 

We search for a general signature of high-mass opposite sign efi pairs and are sensitive 
to both resonant and non-resonant production mechanisms. The CDF collaboration inves- 
tigated the high-mass e/i channel in a Run I search for the direct RPV production and 
decay of the tau sneutrino (ur) . That analysis excluded sneutrino masses below ~ 360 
GeV/c^ for a particular set of RPV coupling values. This Letter describes our continued 
investigation of this channel using data from Run II of the Tevatron taken with the upgraded 
Collider Detector at Fermilab (CDF II). The data sample we use represents an integrated 
luminosity of 344 ± 21 pb~^, approximately three times that used in Run I. We search for 
an efi resonance by examining invariant masses (Mg^) between 100 and 800 GeV/c^ for an 
excess of events. We find event counts that are consistent with SM predictions and use our 
result to constrain models of LFV z/ and Z' decay by excluding e/i coupling values as func- 
tions of the new particle masses. In contrast to searches that utilize low-energy data, where 
the effects of new heavy particles may be masked by cancellations of various contributions, 
the interpretation of our results is robust. 

The CDF II detector is an azimuthally and forward-backward symmetric apparatus de- 
signed to study pp reactions at the Tevatron. The detector has a charged particle tracking 
system immersed in a 1.4 T magnetic field aligned coaxially with the pp beams. A 3.1 m 

n 

long open-cell drift chamber, the Central Outer Tracker (COT) [5J, covers the radial range 
from 40 to 137 cm and provides coverage for the pseudo-rapidity range \ri\ ^ 1 0]. Seg- 
mented electromagnetic and hadronic sampling calorimeters surround thejracking system 
and measure the energies of interacting particles in the range |r/| < 3.6 7, 



A set of drift 



chambers located outside the central hadron calorimeters and another set behind a 60 cm 
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thick iron shield track muons with \ri\ < 0.6. Q|. Additional drift chambers and scintillation 
counters detect muons in the region 0.6 < |?7| < 1.0. Gas Cerenkov counters located in the 
3.7 < 1 77 1 < 4.7 region measure the average number of inelastic pp collisions per bunch 
crossing and thereby determine the beam luminosity. 

We use data collected with high-Py triggers that require central (|?7| ^ 1.0) lepton can- 
didates. We select events that contain at least one reconstructed electron of > 20 GeV 
and one reconstructed muon with > 20 GeV/c. Energy deposited by the leptons in the 
central electromagnetic calorimeters must be of an isolated nature and should also satisfy 
the standard set of "tight" CDF lepton identification (ID) criteria jl^ . 

We select oppositely charged e/i pairs and ensure that the lepton tracks come from a com- 
mon pp interaction. We reject events containing cosmic-ray muons and photon-conversion 



121. In order to maximize 



electrons by imposing additional event topology requirements 
acceptance to a broad range of new physics signatures, we do not impose requirements on 
missing transverse energy or jets. 

SM efi backgrounds include qq —>■ Z/'j* —>■ tt (Drell-Yan), top (tt) and diboson (lyiV, 
WZ, and ZZ) production. We determine geometric and kinematic acceptances for these 



14| based simula- 



processes using the PYTHIA Monte Carlo generator and a GEANT3 
tion of the CDF II detector. These simulations employ the CTEQ5L j3| parton distribution 
functions (PDF's) to model the momentum distribution of the initial state partons. We de- 
fine the acceptance for each background process, a^, as the fraction of generated events 
that satisfy the lepton ID and event topology requirements. We correct the by multiply- 
ing with trigger efficiencies {etrg) measured from W ^ ev and Z —>■ fifi data and factors 
ifreco, fio) that account for differences in lepton reconstruction and ID efficiency between 



simulation and data 



1^. We refer to the combined correction factor, e^^g x 

freco ^ fiDi ^S 6 . 



The expected contribution of each SM background is given by the product of ttj x e 
with the corresponding cross section and the integrated luminosity of our data sample. We 
estimate top and diboson background contributions using next-to-leading order fNLO) cross 
sections: 6.1 pb for ti Q, 12.1 pb for WW, 3.7 pb for WZ and 1.4 pb for |l2|. We use a 
next-to-next-to-leading order continuum {Mu > 30 GeV/c^) cross section for the Drell-Yan 
process, 337.7 pb, which we calculate by scaling the PYTHIA leading order cross section by 
the ratio of NNLO to LO predictions obtained from PHOZPR calculations [iS^ . 

Jets that are misidentified as leptons account for approximately 5% of the total back- 
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ground. We determine the misidentification probability ("fake rate") IJ] by examining 
separate data samples collected with various jet Et triggers for the occurrence of leptons. 
The real leptonic content of these samples is assumed to be negligible. We obtain fake rates 
from the ratios of misidentified leptons and the number of jets in these samples. We then 
estimate the background from this source by applying the measured fake rates, parameter- 
ized as functions of Et, to the jets in events in the high-P^ sample that contain a single 
lepton candidate. 



Channel 


50 < < 100 GeV/c2 


> 100 GeV/c2 




38.8 ± 2.9 


0.6 ± 0.0 


ww,wz,zz 


6.6 ± 0.5 


3.5 ± 0.2 


tt 


3.6 ± 0.5 


3.2 ± 0.5 


Fake Lepton 


2.9 ± 1.7 


^■^ -0.4 


Prediction 


51.9 ± 3.4 


7.7 ± 0.8 


Observation 


56 


5 



TABLE I: Expected and observed event totals for low and high-mass Me/i regions. The uncertainties 
shown are the combination of statistical and systematic errors. 



We present the number of observed and predicted background events for two Me^ re- 
gions in Table HI The dominant uncertainty on the background predictions arises from a 
6% uncertainty in the luminosity measurement Q]. Additional uncertainty contributions 
include those associated with the SM cross sections, fake rate measurements, lepton ID and 
reconstruction scale factors, and PDF model. Overall, we find that uncertainties on our 
signal acceptance and background expectations have little impact on the limits we set. 

The 50 < Mefj, < 100 GeV/c^ region listed in Table H] represents an invariant mass 
range that is rich in background. Finding good agreement between our observation and 
predicted background in this region, we next consider the Mg^ range above 100 GeV/c^. 
Here, too, we find good agreement. Figure ^ shows the observed and predicted background 
Mg^ distributions over a portion of the full 50 — 800 GeV/c^ invariant mass range. 

We quantify the consistency between the M^^ distributions presented in Fig. [T] by per- 
forming a test, using variable- width Me^ bins to achieve sufficient predicted background 
occupancies for the Gaussian approximation of Poisson statistics. The test results in a to- 
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FIG. 1: Observed and predicted Mg^ Distributions. The observed efi invariant mass spectrum 
agrees well with that of the combined SM and fake lepton backgrounds. No events are observed in 
data beyond 159 GeV/c^. 



tal reduced statistic {x^/Ndof) of 14.0/11 = 1.27, under the assumption that systematic 
uncertainties in each bin are completely correlated. The statistic implies a p-value, the prob- 
ability of finding a larger x^, of ~ 23%. This value provides a statistical basis for accepting 
our results as consistent with SM expectations. 

Finding no evidence of new physics in the Mg^ spectrum, we turn from the model- 
independent search to constraints on specific models; RPV sneutrino and LFV Z' decay. 
In RPV supersymmetric models, the s-channel dd — ^{{vi) efi process is governed by 
two LFV couplings. X'^n determines the sneutrino production cross section from d and d 
while Xii2 gives the sneutrino branching ratio to efi |20]. The index i refers to the lepton 
generation of the sneutrino. We do not consider initial-states that include up-type quarks 
since RPV sneutrino hadro-production occurs only through a LiQjDk term in the super- 
potential [2^. The final-state for this process consists of e/x only, i.e., it does not contain 
additional neutrinos or a SUSY LSP. 

Since strong limits exist for z/g and 9^ couplings [2l|, we focus instead on the third 



generation sneutrino, u^. We assume the "single coupling dominance" hypothesis [2lJ and 
set all A and A' couplings but X'^-^^^ and A132 to zero so that contributions to the e/i channel 
originate from the 9^ only. Previous limits on X'^n and A132 from low energy experiments 
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are 0.10 and 0.05 



22j |. assuming squark and slepton masses of 100 GeV/c^. 



We show an example NLO a xBR [2^ that corresponds to X'^^ = 0.10 and A132 = 0.05 
in Fig. 121 We use such curves for different A132 and X'^n values and a mass-dependent axe 
calculated from PYTHIA Monte Carlo simulations of the dd — > z^tI^'t) ~^ e/i process to 
obtain signal expectations over the 100 to 800 GeV/c^ range. Assuming coupling values of 
A132 = 0.05 and X'^^ = 0.1, for example, we find that a hypothetical 300 GeV/c^ Ur signal 
consists of 16.2 ± 1.3 events. 

We calculate an upper limit a xBR for dd — > ^'rl^'r) — >e/i using the numbers of observed 
and predicted background events in bins of Mg^. We scan the Mg^ range in steps of 10 
GeV/c^ and use the simulated mass resolution at each step to weight all data and background 

n 

events between 100 and 800 GeV/c^. A Bayesian algorithm j25| is used with a uniform prior 
to translate the weighted event totals to a 95% confidence level (CL) upper limit on the 
number of signal events in data. We divide this limit by the integrated luminosity and 
signal a X e to obtain a 95% CL upper limit on the a x BR for the process. The axe for 
a generic spin-0 particle decaying to cfi is mass-dependent, increasing slowly from ~ 10% 
at 100 GeV/c^ to ~ 27% at 800 GeV/c^. We account for uncertainties on the signal and 
background expectations in limit calculation although the final results are largely insensitive 
to these numbers. 

The intersection of the observed upper limit and a NLO a xBR curves defines a 9^ mass 
limit for specific values of the X'^^ and A132 couplings, as shown in Fig. |2l We construct 
exclusion regions in the X'^n-Men plane by plotting the mass limit as a function of both RPV 
couplings. Figure El shows the exclusion region parameterized by five assumed values of A132. 
The plot indicates, for example, that we exclude at 95% CL X'^n values above ~ 0.01 for a 
sample i>r mass of 300 GeV/c^ and A132 ^ 0.02. 

Our search is also sensitive to LFV Z' decay. The pp ^Z' ^ cfi process proceeds through 
diagonal f/(l)' couplings at the initial-state vertex and an off-diagonal LFV f/(l)' coupling, 
Q[2, that determines the Z' branching ratio (BR) to efi P]- We set 95% CL limits on as 
a function of the Z' mass using NLO a xBR's from a group of -Ee-inspired models. Although 
Eq models do not incorporate the LFV (^12 coupling by construction, we use these models 
because they provide a theoretical reference by specifying initial-state Z'-quark coupling and 
a NLO Z' production cross section. We utilize NLO cross sections provided by the y, -0 , 
Secluded and rj Eq models (2^ and extend these models to include the (5i2 coupling j22[, 
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FIG. 2: Observed 95% CL upper limit on cr x BR for dd t'T(z^r) e/i (solid line) and the 
NLO prediction (dashed line) as a function of e/i invariant mass. Their intersection gives a 530 
GeV/c^ Dr mass limit for the values of X'^n and A132 indicated. Due to small differences in the 
signal acceptances, the 95% CL upper limit on the Z' a xBR (not shown) is larger than that of the 
sneutrino, ~ 0.02 pb greater at 200 GeV/c^ and ~ 0.003 pb greater at 600 GeV/c^, for example. 



which introduces Z' ^efi decays. 

We set limits on as a function of the Z' mass following the procedure previously 
described for RPV Ur decay. The a x e we use in calculating the upper limit on a x BR is 
obtained from PYTHIA Monte Carlo simulations in which the Z' is required to couple to 
a left-handed leptonic LFV current, as may be favored in i^g-motivated U{iy models that 
incorporate LFV j27|. Z' acceptance (~ 8% at 100 GeV/c^, and increasing to ~ 20% at 
800 GeV/c^) is smaller than that for the scalar sneutrino due to the spin-1 nature of the 
particle. 

The Q[2-Mefj, regions that we exclude for the modified x? 4'^ Secluded and t] models are 
shown in Fig.^J Assuming initial-state Z' couplings specified by the t] model, for example, we 
exclude values above ~ 0.01 — 0.1 for Z' masses between 200 and 700 GeV/c^. Because 
the LFV (5i2 coupling is not intrinsic to Eq models, our limits should not be interpreted as 
constraints on the models themselves. 

In summary, we searched in 344 pb^^ of CDF II data for high- mass e/i events and found 
an distribution consistent with SM predictions. With no evidence for new physics, we 
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FIG. 3: X'^ii — Exclusion Regions. Regions to the left of the five curves shown represent 
ranges of A^j^^^ values that we exclude at 95% CL as a function of Dr mass. Each region corresponds 
to a fixed value of A132. 
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FIG. 4: ~ ^efi Exclusion Regions depicting 95% CL upper limits on the LEV coupling 
in extended ^, ?? and Secluded Eq models. 



set correlated 95% CL limits on the mass and X'^^ and A132 RPV couplings of the z/,-. We 
achieve a ~ 170 GeV/c^ improvement in the mass limit for specific RPV coupling values 
of X'^ii = 0.1 and A132 = 0.05 .2^] and, more generally, exclude z/^ masses over a range of 
Ag]^^ and A132 values. We also place limits on potential LFV couplings of the Z' boson 
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as a function of its mass using E'g-like models of C/(l)' symmetry. 
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